Resonance- and chaos-assisted tunneling in mixed regular-chaotic systems 
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We present evidence that nonlinear resonances govern the tunneling process between symmetry- 
related islands of regular motion in mixed regular-chaotic systems. In a similar way as for near- 
integrable tunneling, such resonances induce couplings between regular states within the islands 
and states that are supported by the chaotic sea. On the basis of this mechanism, we derive a 
semiclassical expression for the average tunneling rate, which yields good agreement in comparison 
with the exact quantum tunneling rates calculated for the kicked rotor and the kicked Harper. 
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Despite its genuine quantal character, dynamical tun- 
neling J] is strongly sensitive to details of the under- 
lying classical phase space |2j. A particularly promi- 
nent scenario in this context is "chaos-assisted" tunnel- 
ing which takes place between quantum states 
that are localized on two symmetry-related regular is- 
lands in a mixed regular-chaotic phase space. The pres- 
ence of an appreciable chaotic layer between the islands 
dramatically enhances the associated tunneling rate as 
compared to the integrable case, and induces strong fluc- 
tuations of the rate at variations of external parameters 
0, 0] . This phenomenon is attributed to the influence of 
"chaotic states" that are distributed over the stochastic 
sea. Since such chaotic states typically exhibit an appre- 
ciable overlap with the boundary regions of both islands, 
they may provide efficient "shortcuts" between the two 
regular quasimodes in the islands [3.0-El Indeed, chaos- 
assisted tunneling processes arise in a number of physical 
systems, e.g. in the ionization of resonantly driven hydro- 
gen , in microwave or optical cavities @ , as well as 
in the effective pendulum dynamics describing tunneling 
experiments of cold atoms in optical lattices 
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While the statistical properties of the chaos-assisted 
tunneling rates are well reproduced by a random ma- 
trix description of the chaotic part of the Hamiltonian 
fl2j |. the formulation of a tractable and reliable semi- 
classical theory for the average tunneling rate is still 
an open problem. Promising progress in this direction 
was reported by Shudo and coworkers 0] who obtain 
a good quantitative reproduction of classically forbidden 
propagation processes in mixed systems by incorporat- 
ing complex trajectories into the semiclassical propaga- 
tor Their approach requires, however, the study of highly 
nontrivial structures in complex phase space, and cannot 
be straightforwardly connected to single coupling matrix 
elements between regular and chaotic states. A comple- 
mentary ansatz, based on a Bardeen-type expression for 
the coupling to the chaos, was presented by Podolsky 
and Narimanov [Tij . In comparison with tunneling rates 
from the driven pendulum, good agreement was obtained 
for large and moderate h, whereas significant deviations 
seem to occur deep in the semiclassical regime 



In the present Letter, we shall point out that nonlinear 
resonances between different classical degrees of freedom 
play a crucial role in chaos-assisted tunneling processes. 
Such nonlinear resonances are known to govern tunneling 
between symmetry- related wells in near-integrable sys- 
tems 0, 0, E3 , where they induce transitions to 
highly excited states inside the well and thereby strongly 
enhance the tunneling rate ^t|- We shall argue that the 
same mechanism is also responsible for the semiclassi- 
cal coupling between regular and chaotic states in mixed 
systems, and determines the average tunneling rate in 
chaos- assisted tunneling. A simple semiclassical expres- 
sion derived from this principle shows indeed reasonably 
good agreement with the exact quantum splittings. 

We restrict our study to systems with one degree of 
freedom that evolve under a periodically time-dependent 
Hamiltonian H (p, q, t) = H(p, q,t + r) and are visualized 
by a stroboscopic Poincare section evaluated at t — nr 
(n e Z). We suppose that H possesses a discrete sym- 
metry which, for a suitable choice of a parameter of H, 
leads to a mixed phase space with two symmetric reg- 
ular islands that are separated by a chaotic sea. We 
furthermore assume that each of the symmetric islands 
exhibits a prominent r:s resonance — i.e., where s internal 
oscillations around the island's center take place within 
r periods of the driving — which manifests itself in the 
stroboscopic section as a chain of r sub-islands that are 
embedded in the torus structure of the regular island. 

The motion in the vicinity of the r:s resonance is ap- 
proximately integrated by secular perturbation theory 
1191 . For this purpose, we formally introduce a time- 
independent Hamiltonian Hq (p, q) that approximately 
reproduces the regular motion in the islands, and de- 
note by (/, 9) the action-angle variables describing the 
dynamics within each of the islands. After the canonical 
transformation ^ d — 8 — s/r ■ 27rf/r to the frame 
that corotates with the resonance, and after averaging 
the resulting Hamiltonian over r periods of the external 
driving, which is justified since i? varies slowly with time 
near resonance, we obtain in lowest nonvanishing order 

H eS {1, 6) = iLzIllfL + 2V r:s cos rd (1) 
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as effective integrable Hamiltonian for the dynamics near 
the resonance. Here, I r:s denotes the action variable at 
resonance, l/m r:s parametrizes the variation of the in- 
ternal oscillation frequency with / at resonance, and V r:a 
characterizes the strength of the perturbation. 

Comparing the pendulum-like dynamics of this effec- 
tive Hamiltonian with the actual classical dynamics gen- 
erated by H provides an access to the parameters of H c s 
without explicitly using the functional form of Hq(j>, q). 
To this end, we numerically calculate the monodromy 
matrix M r:s of a stable periodic point of the resonance 
(which involves r iterations of the stroboscopic map) as 
well as the phase space areas S+ s and S~ s that are en- 
closed by the outer and inner separatrices of the reso- 
nance, respectively. Using the fact that the trace of M r:s 
as well as the phase space areas S^: s remain invariant 
under the canonical transformation to (/,$), we infer 
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I r:s — —(Sr; S + S r . s ), (2) 

— (S^. s — S r . s ) , (3) 
2V 1 

— — = -^-arccos(trM r . s /2) (4) 
mr-a r^r 



from the integration of the dynamics generated by H e ff, 
which allows us to determine I r:s , m r:s , and V r:s . 

The implications of the nonlinear resonance for the cor- 
responding quantum system can be directly seen from 
the representation of the quantized version of H e g in 
the eigenbasis of Hq, which consists of "even" and "odd" 
functions with respect to the discrete symmetry of H. In 
the action-angle variable representation, the eigenfunc- 
tions of Hq are, for a fixed parity, essentially given by 
plane waves ipni'd) ~ exp(irv&) as a function of the angle 
variable, where the integer index n denotes the excita- 
tion as counted from the center of the island. The first, 
"kinetic" term of H e g is therefore diagonal in this basis 
with the matrix elements 



E n = [h(n + 1/2) - I r:s ]7(2m r:s ) , 



(5) 



while the "potential" term 2V r:s cosr$ induces couplings 
between ip n and ijj n ± r with the matrix element V r:a . In 
this way, a perturbative chain is created that connects the 
"ground state" tpo of the island to the excited states ipi r 
with integer I. As was shown in Ref. 0|, this coupling 
mechanism generally leads to a strong enhancement of 
the level splitting between the even and the odd ground 
state in the near-integrable regime, since the unperturbed 
tunneling rate of a highly excited state ipir is much larger 
than that of ipo . 

In the mixed regular-chaotic case, the above tridiag- 
onal structure of the effective Hamiltonian becomes in- 
valid beyond a maximum excitation index n c that marks 
the chaos border, i.e. for which 2irh(n c + 1/2) roughly 
equals the size of the island. Basis states ip n with n > n c 



are defined on tori of Hq that are destroyed by the pres- 
ence of other strong resonances, and therefore exhibit on 
average a more or less equally strong coupling to each 
other. In the simplest possible approximation, which ne- 
glects the presence of partial barriers in the chaos Q , the 
"chaotic block" (H n . n f) n ^ n > >nc of the effective Hamilto- 
nian is therefore represented by a random matrix from 
the Gaussian orthogonal ensemble 0, 0] . 

The probability density P(AE) for obtaining the level 
splitting AE between the ground state energies of the two 
symmetry classes can now be calculated by performing 
the random matrix average over the chaotic part of the 
Hamiltonian. As was worked out by Leyvraz and Ullmo 
[l2j |. this leads to a Cauchy distribution 



P(AE) = 
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with a cutoff at AE ~ 2V c s , where N c and A c denote the 
number of chaotic states and their mean level spacing at 
energy Eg, respectively, and V e s represents the effective 
coupling matrix element between the ground state and 
the chaotic block. In the presence of the nonlinear reso- 
nance inside the island, the latter is evaluated by means 
of the tridiagonal structure within the regular part of the 
Hamiltonian: assuming V r:s to be much smaller than the 
intermediate energy differences, we obtain 
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where the energies Ei r are computed from Eq. (JSJ. The 
elimination of intermediate regular states is performed 
up to the first state -0fer that is already located beyond 
the chaos border [i.e., (k — l)r < n c < kr]. 

Since tunneling rates and their parametric variations 
are typically studied in a logarithmic representation, the 
relevant quantity to be calculated from Eq. © and com- 
pared to quantum data is not the mean value of AE 
(which would diverge if the cutoff is not taken into ac- 
count), but rather the average of the logarithm of AE. 
We therefore obtain the "mean" level splitting AE as 



AE 



V 2 



cir 



N C A ( 



exp 



log- 



N r A r AE 



V 2 



N C A C 



(8) 



where (. . .) denotes the average with respect to the proba- 
bility distribution @ . The expression for the mean split- 
ting further simplifies for our case of periodically driven 
systems, where the eigenphases of the time evolution op- 
erator are calculated. Using the fact that the chaotic 
eigenphases are more or less uniformly distributed in the 
interval < (p < 2?t, we obtain A^ C A C = hu; = 2irh/T. 
This results in the mean eigenphasc splitting 
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FIG. 1: Chaos-assisted tunneling in the kicked rotor at (a) 
K — 2 and (b) K = 3. Plotted are, as a function of iV = 
2tt/H, the splittings Ay> between the eigenphases of the ground 
states at the Bloch phases £ = and n. The associated regular 
islands are shown in the lower panels. The step-like curve is 
the semiclassical prediction Atp for the eigenphase splittings, 
taking into account the 10:2 resonance at K = 2 and the 10:3 
resonance at K — 3. 



FIG. 2: Chaos-assisted tunneling in the kicked Harper at 
(a) t = 2 and (b) r = 3. Plotted are in (a) the splittings 
between the eigenphases of the ground state at £ = and it, 
and in (b) the eigenphase splittings between the symmetric 
and the antisymmetric states that are localized on the center 
of the small, bifurcated islands (at fixed £ = 0) . The step-like 
curves represent the semiclassical predictions, based on the 
8:2 resonance at r = 2 and on the 9:1 resonance at r = 3. 



To illustrate our theory, we apply it to one-dimensional 
systems that are subject to time-periodic kicks. Their 
classical Hamiltonian is given by 

oo 

H(p 7 q,t)=T(p) + rS(t-nT)V(q) (10) 

n— — oo 

where T(p) and V{q) denote the kinetic energy and the 
potential associated with the kick, respectively. The time 
evolution can be represented by the map (p, q) i— > (p, q) 
with p = p — T~V'(q) and q = q + rT'(p), which describes 
the stroboscopic Poincare section at times immediately 
before the kick. The corresponding quantum dynamics 
is generated by the unitary operator 

U = exp (-jT(p)\ exp (~jV{q)j (11) 

where p and q denote the momentum and position oper- 
ator, respectively. 

Specifically, we consider the kicked rotor given by 
T(p) — ip 2 and V(q) = Kcosq with r = 1, and the 
kicked Harper given by T(p) = cosp and V(q) = cosq 
[2fj|. Using Bloch's theorem, we restrict our study to 
eigenf unctions of U that are periodic in position. Wc 
furthermore choose h = 2ir/N where N is an even inte- 
ger. This allows us, for both the kicked rotor and the 
kicked Harper, to write the periodic eigenfunctions as 
Bloch functions in momentum — i.e., with ij)(p + 2ir) = 
ip(p) exp(i£) where tp is the Fourier transform of tp. Since 



the subspace of such functions is TV-dimensional for fixed 
£ |2CJ , the eigenphases and eigenvectors of U can be cal- 
culated by diagonalizing finite N x N matrices. 

Quantum tunneling can take place between a regular 
island in the fundamental phase space cell and its period- 
ically shifted counterparts. As a consequence, different 
Bloch phases £ lead to slightly different eigenphases for 
states that are localized on a given torus in the island. 
The spectral quantity that we discuss in the following is 
the difference Atp — |<^ = °) — tp^ =7T ^ \ between the eigen- 
phases of the island's ground state for £ = and £ = ir. 
(£>(£=°J an d (^(5 =7r ) are calculated by diagonalizing U in 
a suitable basis 0| > and by identifying the ground state 
from the localization properties of the eigenstates near 
the center of the island. Multiple precision arithmetics is 
used in order to calculate splittings below Atp ~ 10~ 15 . 

Figs. Ola), DJb), and[2a) show the eigenphase split- 
tings for the kicked rotor and the kicked Harper, respec- 
tively, as a function of N = 2ir/h, calculated for K = 2 
and 3 in Fig. Q as well as for r = 2 in Fig. Ufa). The 
step-like curves show our semiclassical predictions of the 
eigenphase splittings, which are based on prominent res- 
onance chains boldly marked in the corresponding phase 
space. The relevant parameters m r:s , V r:s and I r[s are 
computed from phase space areas and periodic points 
via Eqs. J2HD- From the numerically calculated phase 
space area S covered by the island, we infer the number 
k of intermediate steps that are necessary to couple the 
ground state to the chaos. An artificially sharp decrease 
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of the semiclassical splitting Aip therefore occurs when- 
ever h passes through a value where S = 2irh(kr + 0.5) 
with integer k. 

In spite of the number of simplifications and approxi- 
mations that are involved in the derivation of the semi- 
classical expression for the mean eigenphase splittings, 
we obtain a relatively good agreement between Aip and 
Aip. In particular, the first major plateau in the quantum 
splittings is remarkably well matched by the semiclassi- 
cal curve, which clearly indicates that the coupling to the 
chaotic sea is mediated by the nonlinear resonance there. 
Our method fails to reproduce the quantum splittings in 
the "anticlassical" limit of large h, e.g. for N < 50 in 
Fig-H^b). Preliminary calculations show, however, that 
a better agreement in this regime might be obtained by 
properly taking into account the action dependence of 
Vr; S in the effective Hamiltonian QJ, which is completely 
neglected in the present treatment. More details will be 
presented in a subsequent publication. 

Apart from Fig.QJb), where also plateaus of higher or- 
der are well reproduced by the semiclassical theory, we 
observe a systematic tendency to overestimate the ex- 
act quantum splittings for low and moderate values of 
l/h. We tentatively attribute this fact to the existence 
of partial barriers in the chaotic part of the phase space, 
which may enhance the effective size of the island for the 
quantum tunneling process. In particular, it is known 
that "Cantori" (i.e., broken tori) in the chaos inhibit the 
quantum transport in a similar way as invariant tori in 
the island, as long as the phase space area associated with 
the classical flux through the Cantorus is smaller than 
ttH \'2\\ . A better agreement with the quantum splittings 
might therefore be obtained by properly incorporating 
hierarchical states |22| , which are localized in the imme- 
diate vicinity of the island, into the semiclassical descrip- 
tion (see in this context also @). 

Finally, Fig. 0b) shows the case of tunneling between 
two symmetric regular islands in the kicked Harper at 
t = 3, which arise from a bifurcation of the central island 
taking place at r = 2. The quantum splittings are now 
given by the eigenphase difference between the symmetric 
and the antisymmetric state associated with the pair of 
islands, calculated here at fixed £ = 0. We see that the 
splittings display a prominent plateau at N ~ 300 . . . 500, 
which is well reproduced by the semiclassical prediction 
based on a 9:1 resonance inside the islands. 

In conclusion, we have presented a straightforward 
semiclassical scheme to reproduce tunneling rates be- 
tween regular islands in mixed systems. Our approach 
is based on the existence of a prominent nonlinear reso- 
nance inside the island, and uses elementary classical pa- 
rameters associated with this resonance to estimate the 
coupling rate from the island to the chaos. In combina- 
tion with a random matrix description of the chaotic part 
of the Hamiltonian, we obtain a simple expression for 
the average level splittings between symmetry-related is- 



lands, which agrees reasonably well with the exact quan- 
tum splittings calculated for the kicked rotor and the 
kicked Harper. Our study underlines that nonlinear res- 
onances govern the coupling between regular islands and 
the surrounding chaotic sea in the semiclassical limit. We 
expect that they play an equally prominent role also in 
complex systems with more degrees of freedom, for which 
our approach could develop into a useful method to quan- 
titatively estimate tunneling rates in presence of chaos. 
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